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ABSTRACT 


In this research paper, the Laplace transform method combined 
with the semi-analytical Adomian decomposition method 
(LADM) is proposed to solve the mathematical model of crime 
deterrence in society. The model is solved to obtain analytical 
solution to the governing parameters in the form of a rapidly 
convergent series to illustrate its reliability, capability, and 
efficiency of this hybrid method. The practical result obtained 
reveal, the method is accurate and an efficient tool for solving 
wide variety of several first and higher order models. 
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INTRODUCTION 


The term global warming is defined as the unabated 
steady rise in the level of atmospheric carbon (iv) 
oxide gas beyond an allowable threshold. From pre- 
industrial times, the level of CO2 was just 270ppm 
(Parts per million). This number has risen to a 
worrisome level of 395ppm_ (Parts per 
million).Misra and Verma (2013), Rasool and 
Schneider (1971), (Tennakone, 1990). This rising 
trend of atmospheric carbon (iv) oxide if not 
checked does not bode well for both human and the 
ecosystem balance. Sequel to this, group of 
industrialized nation under the aegis of G7 have 


gathered severally to brainstormed on how to 
mitigate the effect of climate change occasioned by 
the ominous rise of carbon (iv) oxide (Pluss, 1956), 
McMichael, Woodruff and Hales (2006). Concerted 
efforts have been made by these countries to reduce 
and stabilize the future concentrate C0, via 
legislation, Woodwell, Hobbie, Houghton, Melillo, 
Moore, Peterson, and Shaver (1983). Reports 
abound that shows climate change accounts for 
about 90% of health risks or challenges faced by 
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human. Therefore, if this rise in carbon (iv) oxide is 
not checked, it will portend danger to humans. For 
example, extreme weather conditions such as 
floods, Tsunami, windstorm, and drought are all 
caused by climate change. The attendant effects of 
the above conditions are direct injuries, 
malnutrition, infectious diseases, and airborne 
disease, Marten, Jettsen, Niessen and Rothmans 
(1995). 


Equally, increased in the presence of vectors during 
this period cause vector-borne disease such as 
malaria, diarrhoea, dengue fever and _ others 
(Alexiadis, 2007). The accompanying heat waves 
also evoke respiratory and cardiovascular problems 
to humans. Following the above adverse effects of 
carbon (iv) oxide in the atmosphere, a better 
understanding behind the factors responsible for the 
increased level of CO, in the atmosphere and their 
effects is needed. Previous studies have shown that 
forest biomass and human population are the 
primary factors responsible for the rise in the level 
of carbon (iv) oxide. Human activities like 
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deforestation and burning of fossil fuel triggers the 
level of atmospheric C0. Upstaging the forest- 
biomass carbon (iv) oxide link due to human 
activities also lead to increase in the level of 
atmospheric carbon (iv) oxide, Detwiler and Hall 
(1988), Khasnis and Nettleman (2005), (Kurane, 
2010). 


Volume of studies have been devoted to 
understanding the interaction of human population, 
carbon (iv) oxide and its attendant effects on human 
population and climate change. Caetano, Gherardi 
and Yoneyama (2011) have considered the 
mathematical model that relates atmospheric 
C0,, forest biomass and GDP. In this study, it was 
shown that clean technology and reforestation are 
needed to attain desired C02 level. Malhi and Grace 
(2000), (Ewers, 2006), Kremen, Niles, Dalton, 
Daily, Ehrlich, Fay, Grewal and Culley (2000), 
(Mahar, 1989), Olabisi, Reich, Johnson, 
Kapuscinski, Suh and Wilson (2009) have studied 
the connection between human activities and global 
warming using feedback control. The study reveal 
increased in the level of carbon (iv) oxide has 
destabilizing effect. The biomass-carbon (iv) oxide 
system have been investigated using a mathematical 
model. The study show that excessive deforestation 
destabilizes the biomass-carbon (iv) oxide 
equilibrium. 


The Laplace transformation and Adomian 
decomposition method (LADM) is a hybrid semi- 
analytical method which is a fusion of the two 
methods. The method requires taking the Laplace 
transform of both sides of the equation using the 
accompanying initial condition. Thereafter, the 
resulting solution is then written in operator form 
which gives an n —fold integral proportional to the 
highest degree of the invertible function. The zeroth 
order and the recursive algorithm are then obtained 
which gives the solution that converges to the exact 
solution if it exists. The LADM originally proposed 


by Khuri, Songen and Mendelson (2003), (Khuri, 
2004), Khuri and Alchikh (2019, 2020) is 
preferable over the famous Adomian decomposition 
method because it converges faster to the exact 
solution. It has been extensively applied in the 
following areas: linear and nonlinear PDEs, 
coupled systems of PDEs, numerical solution of the 
Duffing equation, analytical solution of an HIV 
model, Newell-Whitehead-Segel equation, systems 
of ordinary differential equations, linear and 
nonlinear integral equation with weak kernel, nth- 
order integro-differential | equations, = two- 
dimensional viscous fluid with shrinking sheet, 
logistics equation, numerical solution of the crime 
deterrence model, nonlinear function equation and 
convection diffusion-dissipative equations, (Fadaei, 
2011), Khan, Hussain, Jafari and Khan (2010), 
(Yusufoglu, 2006), (Nasser, 1997), (Ongun, 2011), 
(Pue-on, 2013), (Cherrualt, 2002), (Dogan, 2012), 
(Hendi, 2011), (Wazwaz, 1999, 2010), (Waleed, 
2013), (Manafianheris, 2012), Mohamed and Torky 
(2013), Koroma, Widattala, Kamera and Zhang 
(2013), Islam, Khan, Faraz and Austin (2010), 
Yindoula, Youssouf, Bissanga, Bassino and Some 
(2014), (Khuri, 2001), Al-Khaled and Allan (2005), 
(Doan, 2012). 


In this present study, we seek analytical solution to 
the parameters governing the problem using 
Laplace Adomian decomposition method. The 
article is composed as follows: section 1 takes the 
introduction. In section 2, the basics, and 
fundamentals of the Adomian decomposition 
method is extensively discussed. Section 3presents 
an in-depth review of the hybrid Laplace Adomian 
decomposition method. In sections 4 & 5, the Padé 
approximation and application of LADM to the 
model problem is presented. Numerical application 
via simulation is contained in section 6 and finally 
the section 7 gives the concluding remarks. 
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ADOMIAN DECOMPOSITION METHOD (ADM) 
Consider a general nonlinear differential equation of the form 


Fly@)] = g@) (1) 
Where F is a nonlinear operator and y, g are both functions of x 

Decomposing the nonlinear operator, F into two parts comprising a linear and nonlinear operator. 
Lly@)] + RIy@)] + N[y@)] = g@) (2) 


Where L is the highest order derivative that’s assumed to be invertible, R is the linear differential operator 
with order less than that of L, N is a nonlinear term and g is the source term 


Rewriting Eq. (2) for L[y(x)], we obtain 


Lly(x)] = g@) — RIly@)] — Niy@)] (3) 
Taking the inverse operator, L~* on both sides of Eq. (3), we get 
y(x) = L*g(x) — LR y(@)] - LN [y(x)] (4) 


Where ¢ is the term arising from the integration of the source term.It is obtained using the following 
sequence depending on the order of the given equation. 


d 
y(0) forL= aT 
2 
y(0) + xy'(0) for L == 
ba 3 
_ y(0) + xy'(0) + zy 0) fork. = ax 
= x2 i d* 
y(0) + xy (0) + rd (0) + 31Y (0) for L = aa 
! a " x" hn pt (n) dn** 
y(0) + xy OO) Se O+Ty CO) tarry (0) jor == ag 


By the standard Adomian decomposition method, we write the unknown solution as an infinite 
decomposition series of the form 


y(x) = Un=0 Yn) (5) 
Putting Eq. (5) into Eq. (4), we obtain 
Ln=0 Yn(X) = & — LT RD n= Yn] — LN Din=0 Yn] (6) 
Matching both sides of Eq. (6), we obtain the zeroth order component given by 
yo = 
Then the recursive relation is given by 
Yna1(*) = —L*R[yy] ~ L“*N[yn],n 20 (7) 
The solution of the problem in Eq. (1) is obtain as limit of the decomposing series 
y(x) = limp + Yn (x) (8) 
Similarly, the nonlinear term can be determined by an infinite series of the Adomian polynomials. That is, 
NlY¥0.¥1» Yar ++ Yn] = Ln=0 An (9) 
Then the A; are obtained from the relation 
a” ee 

n= Sam IN Dio Ayaan = 0,1,2,3 (10) 
Using Eq. (9), the first five Adomian polynomials are given as 
Ao = N(¥) 
A; = YiN' (yo) 
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Az = Y2N'(¥o) + avi N" Oo) 
1 
Az = y3N'(¥o) + V1v2N" (Yo) + zi No) 
’ 1 ”" 1 mw 1 (iv) 
Ay = ¥4N'(¥o) + aN (yo) (213 + ¥3) +5 aN Wo)¥iv2+— A N“” (yo) yt 


if 1 4 
As = ysN'(yo) + 5N (Yo) (2y1¥4 + 2V23) +z iN (yo) Byfy3 + 3y1y5) +7 aN ©) (yy) (ye y¥2) 


1 Vw) 
ta (Yo) Yi 


1 1 
Ag = YeN' (yo) + aN" On) (219s + 2y1y4 + y3) +3 aN" Oo) Byr V4 + V3 + 6Y1¥2Y3) 


1 (iv) 3 5 (v) 1 (vi) 6 
+ NO Oo 4vi ya + Oi ¥2) + NY oye + NOY Oo) 1 


1 it 
Az = y7N'(yo) +5 SN" W0) (291%. + 2V2¥s + 2y3Va) Lay ETL (yo) (3yi ys + 3y1¥s + 3By3ys + OV V2V4) 
- 1 
ay (yy) (4y?y4 + 12y? yoy3 + bas ig (yo) Sytys + 10yzy3) 


1. 1... 
+78 Ci) (y))yPy. + he iD (y))y] 


1 
Ag = YeN'(Yo) +N" Wo) 2y1¥7 + 2¥2Vo + 2¥a¥s + ye) 


1 
tae (yo) (3y2 V6 + 3¥aVE + 3¥2V3Z + 6V1V2Vs + 6¥1¥3V4) 


— 
+ FN Oo Ayiys + 12yi Vays + Lyi viva + Ovi ys +2) 

1 1 

+ TN (Vo) Syiva + 20y?y2¥a + LOyZy3) + FN? Wo) Oris + 15yty3) 


7 1 
ti Mt) (y))y Sy. + ad (tid (y,) y8 


LAPLACE ADOMIAN DECOMPOSITION METHOD (LADM) 
In this subsection, we outline the basics of the fundamentals of the fusionLaplace transformation and 
Adomian decomposition method (LADM) 


Consider a functional differential equation of the form 


L[u(x)] + R[w@x)] + N[u@)] = g@) (11) 
Subject to the initial condition 

u(x,0) = f(x), “2? = ace) (12) 
Rearranging the ae we obtain the following relation for L[u(x)] 

L[u(x)] = g(x) — R[u@)] — N[u@)] (13) 


Applying Laplace transform on both sides of Eq. (11), supposing the highest differential operator is of order 
two and using the differentiation property, we get 


s*L{u(x)} — sh(x) — f(x) = L{g(x)} — L{Ru(x)} — L{Nu(w)} 
s*L{u(x)} = sh(x) + f(x) + L{g(x)} — L{Ru(x)} — L{Nu(x)} 


L{u(x)} = "2+ 24S cf{g(a)} - SL{RuG)} -— 5 L{Nu(x)} (14) 
Next, we ee the inverse transform on both sides of Eq. (14), we obtain 
u(x) = $x) ~ £74 | L{Ru(x)} — 5 L{Nu(x)}| (15) 


Where #(x) is the term arising from the first three terms on the right-hand side of Eq. 
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Next, we assume the solution as decomposing series in the form 

U(X) = Yin=o Un(X) (16) 
Similarly, the nonlinear terms are written in terms of the Adomian polynomials 

Nu(x) = Yn=04n (17) 
Where the A;? represents the Adomian polynomials defined in the form 


= > INS A'N) Ing = O23 (18) 


n ~ nldaan 


Plugging Eqs (16) and (17) into Eq. (15), we obtain 


Dito Un (x) = G(x) — £72 |S LER Dio ta (x)} -— GLEN Dieveo An} | (19) 


Matching both sides of Eq. (19), we obtain an iterative algorithm in the form 


Up (x) = h(x) 
1 . 1 . 


up (x) = —£72| 5 LER Dito ta (2)} — GLEN Davao As} | (20) 


, ‘i : sco} — =f I 3 a 


u,(x) = -—L7} 


ie(xy= = * 


Un+1 (x) = —£ 


nya 


n=0 


1 = 1 
zee > nc ~ Peng 


n=0 


Then the solution of the differential equation is obtained as the sum of decomposed series in the form 
U(x) © Up (x) + u(x) + uz2(x) + (21) 


PADE APPROXIMATION 

In Mathematics and other applied sciences, power series representation of a function is usually in truncated 
form. To approximate these functions to an appreciable degree, polynomials are used because their 
singularities are easily noticeable in each finite region. However, the radius of convergence may not be large 
enough to contain two boundaries, for this reason, power series is not always the best method to 
approximate a function. To overcome this inherent hurdle, a new approximation is applied to the solution 
obtained using power series as a quotient of two functions with varying degrees in a finite interval. 


Padé approximation has been widely used to approximate several problems and has tremendous applications 
especially in computer calculation because it gives a better approximation without truncating its power 
series and still in problems where the series diverges. The different Padé approximants are obtained with the 
use of symbolic software Mathematica 


A rational approximation to a function f(x) on [a, b] is the quotient of two polynomials, Py(x) and Qy(x) 
of degrees N and M respectively. It is denoted by [V/M](x) 


That is, [V/M](x) = a ,a<x<b (22) 
M 

Now consider the formal power series 

f (x) = Dito Cex" (23) 
_ Pn) N+M+4+1 

f(x) = 282 + O(N M44) 

Rearranging gives 
_ Pu) _ N+M+41 

f(x) — SBD = O(eNMH) (24) 
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Multiply both sides of Eq. (24) by a constant keep it unchanged, hence we impose the normalization 
condition. 

Qu (0) = 1.0 (25) 

Next, we require that Py (x) and Qy(x) have non-common factors, so we write the coefficients of the Py (x) 
and Qy(x) as follows 

Py(X) = Ap + a,x + Ayx? +2 + ayx% 


26 
Qu (x) = bo + byx + box? tet by x! ( ) 


Using Egs. (25) and (26), we multiply Eq. (23) by Q(x). This linearizes the coefficient equation. It is given 
in detailed form as 
(Ap + a4X + Agx? +++ ayx) = (bo + yx + box? +2 + dy x!) (co + yx + Cox? + + Cy x) 


N+1 ,N+2 
, yee 


Equating the coefficients of x“*1, x ., xN*™ successively to zero, we obtain the system of Equations 


byCn—m+1 + Dy-1eN-m42 + °° + Doty 41 = 0 
buCn-m+2 + by-1Cn-m+3 + + + Docn+2 = 0 (27) 
byCy + by-1¢n41 + °** + DoCn+m = 0 
For j < 0, we define c, = 0 for consistency. Setting by = 1, Eq. (22) become a set of Mlinear equations for 
M unknown coefficients in the denominator. 


Cn-M+1 Cw—Mm+2 Cnt 1 be Cy 
| = |e 


Cu+m 


a 


Cn-m+2 ©N—M43.......C v4.2 


N 


bm-1 


mie | | (28) 
Cy ONG wicseiniicanccc NG MGA by 


Solving the above system in Eq. (28), the coefficients b; for i= 1,2,....M may be found. Since the 
coefficients of the numerator, Co, C1, C2, ...., Cx is known. We equate the coefficients of 1, x,x?,...,x%,xNt™ 
to obtain the remaining coefficients A, A1, Az, ..., Ay 

x° sc, ofdyend 

x1 ucyb, +ieesearch O 

x: Cob2 + Cb, + C2 — AZ = 0 

XP B.Co bay PUD GP ely — ay G0 

x: Coby + Cy by_-1 + Cy — An = 0 


Writing the above in explicit form, we obtain 


Ap = Co 

a, = C1 + byCo 

Ay = Cy + b,c, + boc, (29) 
po 2°0 


Ay = Cy + Lees DeCw-n,N = min(N, M) 


In view of Eqs. (28) and (29), the numerator and denominator of the Pade approximant are all determined 
which agrees with the original series to the order of x“*™- 


Now, to solve the system (29) for the set of unknowns, we assume that the Eqs. (28) and (29) are non- 
singular, so we can obtain the solution via the determinant. See [42-45] 


an-M+1 QN-M+20 ee ant+1 
an On+1 PRETEE an+M 
NW od ht SNC cae ay eres 
ee eae ee (30) 
~~ QN-M+1 &N-M+2---- an+1 
an ANF iseserssesee aN+M 
aM M-1 1 


To obtain the diagonal Pade approximants of different orders such as [2/2], [4/4], [6/6], we use symbolic 
software Mathematica 


APPLICATION OF LADM TO THE MODEL 
Given the mathematical model for global warming as follows 
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= Qo + AN — aX — A, XF G1) 
aN N 

@ = sN(1-~)-6XN +nNF (32) 
aF F 

<= UF (1 = =) — NF +1,A,XF (33) 
Rearranging the above Eqs(31) — (33), we obtain 

= Qo + AN — aX — A,XF (34) 
= = =N(L—N) —OXN + ngNF (35) 
<=" F(M — F) — ONF + 17,A,XF (36) 
X(0) > 0,N(0) > 0,F(0) >0 (37) 


Where the parameters X, N and F have their usual meanings 
X = amount of C02 in the atmosphere 

N = Available human population 

F = Forest Biomass 


Taking the Laplace Transform of both sides, we get 


e {=} = L{Qo} + L{AN} — L{aX} — LIA, XF} (38) 
Z {=} = = £{N(L — N)} — L{OXN} + L{npNF} (39) 
7 {=} = £{F(M — F)} — L{@NF} + L(A, XF} (40) 
Applying the formula for Laplace transform in the first derivative, we obtain 

wLl{X} — X(0) = L{Qo} + AL{N} — al{X} — A, L{XF} (41) 
wL{N} — N(0) == L{N(L — N)} — OL{XN} + WPL{INF} (42) 


wL{F} — F(0) = -L{F(M — F)} — PL{NF} + 1A, L{XF} (43) 
Using the initial condition in Eq. (37), Eqs (41) - (43) reduced to 


WL{X} = oe + AL{N} — aL{X} — AL{XF} (44) 

wL{N} = > £{N(L — N)} — OL{XN} + nbL{NF} (45) 

wL{F} = + L{F(M — F)} — PL{NF} + 1,4, L{XF} (46) 

Further simplification by dividing both sides by w, we get 

L{X} = 224 cn} - <c{x}- 2 £fA} (47) 

L{N} = = LEN(L — N)}-—£{B} + cfc} (48) 

L{F} =~ c{F(M — Fy}-2 cfc} 4+ £f4} (49) 

Where A = XF,B = XN,C = NF (50) 

By the Laplace transform decomposition method, we represent the solution as infinite series of the form 
X = Yn=0XnN = Un=o Nw F = Un=o Fn (51) 


Where the terms X, N,, and F, are to be obtained via the recursive relation. Similarly, the nonlinear 
operators, A,B and C are decomposed as follows 

A = Yn=04n B= Yn=0 Bu C = Yn=0 ln (52) 

Where A, B, and C,, are the Adomian polynomials. The first eight of these polynomials are given by 

Ag = Xoo 

Ay = XoF, + X1Fo 

An = XoF, + X1F, + X2Fo 

A3 = XoF3 + X1F2 + X2F, + X3F 
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Ag = XoFy + X1F3 + X2Fo + X3F, + X4Fo (53) 
As = XoFs + X Fy, + X2F3 + X3Fo + X4F, + X5Fo 
Ag = XoFg + X1F5 + XoF, + X3F3 + X4Fo + X5F, + Xe6Fo 
Az = XoF, +X Fe + XoF5 + X3F, + X4F3 + X5Fo + X6F, + X7Fo 
Ag = XoFg + XyF7 + X2Fo + X3F5 + X4Fy + X5F3 + X6Fo + X7F, + XgFo 
Bo = XoNo 
By, = X9pN, + X1No 
By = Xo Nz + XN, + X2No 
Bz = X)N3 + X1,Nz + X2N, + X3No 
By = XoNq + X1N3 + X2N2 + X3N, + X4No (54) 
Bs = XoNg + X1Nq + X2N3 + X3N2 + X4N, + X5No 
Be = XoNg + X1N5 + X2Nq + X3N3 + X4,Nz + X5N, + X6No 
Bz = XoN7 + X,Neo + X2N5 + X3Nq + X,N3 + X5N2 + X,6N, + X7No 
Bg = XoNg + X,N, + X2No + X3N5 + X4N, + X5N3 + X6N2 + X7N, + XgNo 
Co = NoFo 
Cy = NoF, + Ni Fo 
Cz = NoF2 + Ni Fy, + N2Fo 
C3 = NoF3 + Ni Fp + NoF, + N3Fo 
Cy, = NoFy + N,F3 + NoF2 + N3F, + N4Fo (55) 
Cs = NoFs + Ni Fy, + NoF3 + N3F, + NaF, + NsFo 
Ce = NoFe + Ni Fs + NoF, + N3F3 + NgF2 + NsF, + NeFo 
Cz = NoF7 + Ni Fo + NoFs + N3F, + NaF3 + Ns F2 + NoF, + N7Fo 
Cg = NoFg + N,F, + NoFo + N3Fs5 + NyFy + NsF3 + NoF2 + N7F, + NeFo 
Putting Eqs. (51) and (52) into Eqs (47) - (49) 


LAD Xn} = SS + LEN Nn} — LER Xn} — SLD0 An} (56) 
co co CO 6 co co 
LAD Nn} = — LAD 9 Nn (L — Dito Nn)} — SLD Bn} + PLO 0 Cn} (57) 
co co co co A co 
LAD 0 Fab = LID Fa (M — Lg Fa} — SLD Cn} + ALD An} (58) 
Matching both sides of Eqs (56) — (58) yield the following iterative algorithms 
Qo 
L{Xo} = we 


£{Xq} = LEN} — <£{X%o} - 2 LfAg} 

£{X,} = = £{N,} — < £644} - 2 {Ay} (59) 
L{Xg} = = L{N,} — = £{X,} - 2 £{AQ} 

£Xpas} = ~ £{Nq} — {L{Xn} — Lf An} 

L{No} = 0 

L{N,} = = L{No(L — No)} — —£{Bo} + “2 £{Co} 


L{N2} = — LEN, (L — Ny)} - = £{B} + £6} (60) 
L{Ng} = = L{N2(L — Nz)} — = £{B2} + 2 £{C,} 
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Ss 6 Tp 
LENn+13 = — etn — Nn) — EBS - — LtCnt 


L{Fo} = 0 
L{F,} = £EFy(M — Foy} - 2 £{ Co} + = £fAg} 
£{F,} =~ LF, (M — F,)}- 2 £{C,} + cay} (61) 
_ Hu p 1A, 
L{F,} = —— L{F,(M — F2)} - — £{C,} + = £fAy} 
Fae = com — 8} - © ef} + hg) 
wM w w 


NUMERICAL APPLICATION 
In this section, we apply the LADM to the numerical solution of the model using simulation. 


Applying the inverse Laplace transform to both sides of Eqs (59) — (61) gives 


Q 1 1 
Substitution of Eq (62) into the second equations in Eqs (59) — (61), we get 
A AQ  AxQo 
AMS a ye to y3 
s(WL- ) QQ Tp 
L{N,} = ea Ww + a, Z, ‘ (63) 
uw (WM — 14410 
8) = Sala) as ts 
Putting the values of L{X,}, L{N,} and L{F,} into the second Eqs. (59) — (61), we obtain 
L _ A AQ A1Qo 
ad = oa ye to y3 
(WL-1) — QQ 
L{Np} = SO — S00 4 7P (64) 
_ BH (WM — i P@  1A1Qo 
UF) = Sg ae) ps ta 


Evaluating the Laplace transform of the quantities on the RHS of Eqs. (62) — (64), and applying the inverse 
Laplace transform, we obtain the values, X,(t), N,(t), F,(t()and X>(t), N2(t), F,(t). Similarly, the other 
higher order solutions X3(t),X4(t),....,X,(t), N3(t), N4(0),...., Ny (0), F3(t), F, (0), ...., F(t) are obtained 
recursively in a similar fashion using Eqs. (59) — (61) 


Now to obtain the solution of the parameters of interest in explicit form, we apply LADM to the model by 
taking the following values via simulation. We take X(0) =1,N(0) =1,F(0) =1, for the three 
components of the model. Next, we take Q) = 1,4 =0.05,a@ = 0.03,2, = 0.0001,s == 0.01,L = 
1000, 8 = 0.00001, v = 0.2,M = 2000, = 0.01, d = 0.0002, 7, = 0.01. A few first approximations for 
X(t), N(t) and F(t) are calculated and presented below using LADM as follows. 


£{Xo} = , L{No} = — ,L{Fo} = — (65) 
L{X,}= ae 7 — = 
WwW WwW WwW 
L{Ny} = — a —— a —_ (66) 
0.00999 0.0002 0.000002 
A ee 


Taking the inverse Laplace transform of both sides of the above equations, we obtain the solutions of the 
parameters as follows. 
X(t) = 0.05t — 0.03t? + 0.0000333333t° 
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N(t) = 0.00999t + 10° °t? — 3.33333 x 10° °t? (67) 
F(t) = 0.00999 — 0.0001t? + 3.33333 x 1077t? 
Next., we calculate the [5/5] Pade approximates of the infinite series solution which gives the following 
rational fraction approximation of the parameters of interest using Mathematica 
1. + 0.7195405t — 0.365002t? + 0.0050¢? — 0.0000051t* 
Xpade(t) = 1 


— 0.593t + 2.16 x 1071 t? +1.68 x 10-15¢3 — 1.4 x 10-16¢4 


1.+0.0049t—0.048t2—0.0333 34+1.6649x1078t* 
1-0.0045t+6.0408x10-17t24+4.74x107-1 ¢34+4.9x10-19t4 


ae 1. + 0.0502t — 0.00506t? + 0.05308t? — 1.6565 x 10-8t*# 
Pade\’” ~~ 1 —0.049t + 6.9697 x 10-142 — 3.503 x 10-17t3 — 3.4 x 10-18¢4 


Npade (t) = 


RESULTS AND DISCUSSION 

In this subsection, the results of the problem in Eq. (1) are presented to show the effects of the governing 
parameters on the model. The effectiveness and accuracy of the numerical methods are displayed in Tables 
1-3 and Figures 1-7. The methods give highly accurate results in few steps. The results obtained when 
compared are consistent with literature 


Table 1: Numerical Computations for X(t) 


FAD OFA D)\Y Gis OVW D)\Y oT: (Se beak Oct (3) oS 
1 


0 1 1 
0.2 | -1.27322 -1.27369 -1.27310 
0.4 | -13.01120 -13.1023 -13.1201 
0.6 | -39.11090 -41.6953 -41.6102 
0.8 | -92.2623 -134.119 -134.102 
1.0 | -196.788 -195.23 -194.21 
-392.44 178.808 178.801 


Table 2: Numerical computations for N(t) 
t LADM LADM-Padé 4" Order R-K 


0 1 1 1 

0.2 | 19.0903 19.0903 19.0901 
0.4 | 37.3871 SZ 37.3862 
0.6 | 56.0688 56.0699 56.0700 
0.8 | 75.5955 75.6116 75.6110 
1.0 | 96.867 96.9949 96.9950 
[2 | 121.397 122.11 122.020 


Table 3: Numerical Computations for F(t) 
LADM LADM-Padé | 4" Order R-K 


0 -1.05 -1.05 -1.05 
0.2 | -1.04715 -1.04715 -1.04715 
0.4 | -1.02824 - 1.02796 - 1.02796 
0.6 | -0.983335 -0.973338 -0.97337 
0.8 | -0.959682 -0.833489 -.0.833452 

: -1.30508 -0.318347 -0.318340 
1.2 -3.3568 -5.94349 -5.94340 
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— LADM 
— LADM-Pade 
— RK4 


x(t) 


Time 


Figure 1 Computation of Atmospheric CO, Against Time 


— LADM 
— LADM-Pade 
— RK4 


Nit) 


-4 ie! -3 —3 -1 a 1 ra 


Time 


— LADM 
— LADM-Pade 
— RK4 


P(t) 


-5 —4 -3 -2 -T 0 4 2 


Time 
Figure 3 Computation of Forest Biomass Against Time 
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= — Atmospheric C02 
~~ Human Population 
-i1.0 -0.5 0.0 0.5 1.0 14 2.0 
Nit) 
Figure 4 Computation of Amount of CO2 Against Human Population 
— Biomass Forest 
= ~~ Human Population 
-5 -4 -3 —2 -1 o 1 rs 
F(t) 
Figure 5 Computation of Amount of CO2 Against Forest Biomass 
0.20 
0.715 
0.70 
= —— Human Population 
> 0.05 
7 ~~ Biomass Forest 
0.00 
-0.05 
-0.10 
-0.5 0.0 05 1.0 14 2.0 
Fit} 


Figure 6 Human population Against Forest Biomass 
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Figure 7 Computation of Atmospheric CO2, Human Population and Forest Biomass against Time 


CONCLUDING REMARKS 

In this study, the approximate analytical solution of 
the mathematical model describing the dynamics of 
carbon dioxide in the atmosphere is solved using 
the fusion of Laplace transform and Adomian 
decomposition method (LADM). The validity, 
accuracy, flexibility, and effectiveness of the 
method is demonstrated by obtaining the exact 
solution of the parameters of interest subject to the 
initial condition. The solution obtained shows the 
MADM is effective and convenient. 


Furthermore, MADM is a promising tool to 
effectively both linear and nonlinear PDEs. The 
benchmark solution is a ready reference for further 
works in the crime model. 
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